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ABSTRACT 

We present deep integrated-light spectroscopy of nine M33 globular clusters taken with the Hec- 
tospec instrument at the MMT Observatory. Based on our spectroscopy and previous deep color- 
magnitude diagrams obtained with IIST/WFPC2, we present evidence for the presence of a genuine 
intermediate-age globular cluster in M33. The analysis of Lick line indices indicates that all globular 
clusters are metal-poor ([Z/H] < —1.0) and that cluster M33-C38 is ~5— 8 Gyr younger than the rest 
of the sample M33 star clusters. We find no evidence for a population of blue horizontal branch stars 
in the CMD of M33-C38, which rules out the possibility of an artificially young spectroscopic age due 
to the presence of hot stars. We infer a total mass of 5 — 9 x 10^ Mq for M33-C38, which implies that 
M33-C38 has survived ~ 2 — 3 times longer than some dynamical evolution model predictions for star 
clusters in M33, although it is not yet clear to which dynamical component of M33 - thin disk, thick 
disk, halo - the cluster is associated. 

Subject headings: globular clusters: general — galaxies: star clusters — galaxies: evolution — galaxies: 
formation — galaxies: structure 



1. INTRODUCTION 

Star clusters which have survived for several (~ 2 — 8) 
billion years (intermediate-age clusters) can provide im- 
portant clues concerning both the evolution history of 
their parent galaxies and for understanding the destruc- 
tion processes which erode star cluster systems (e.g., 
iGoudfrooii et alJ 12004 and references therein). In the 
Galaxy, clusters of this age all fall i n the categ ory of "old 
open clusters" (see compilation in lFriellFlQQlTD - they all 
reside in the Galactic thin disk, and have lower masses 
('^few 10^ Mq) than typical ancient globular clusters 
(> lO^-S Mq). The oldest of these old open clusters help 
to pinpoint the age of the thin disk, while the ages and 
masses of the entire population help constrain the sur- 
vival rates of clusters in the Galactic disk. Intermediate- 
age clusters are also known to exist in the Large and 
Small Magellanic Clouds (LMC and SMC). The LMC 
contains a population of clusters which are 1 — 3 Gyr, 
while the SMC formed numerous clusters 4 — 8 Gyr ago. 

In more distant galaxies, it becomes increasingly dif- 
ficult to establish the presence of intermediate-age clus- 
ters, primarily because the techniques available to study 
compact clusters become more limited. For example, 
the integrated optical colors of intermediate-age clus- 
ters are degenerate with those of ancient (> 8 Gyr) 
clusters. Although a number of clusters have absorp- 
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tion line index strengths (measured from integrated-light 
spectroscopy) which suggest that these objects are of in- 
termediate age, their absolute ages cannot be definitively 
established due to the possible presence of hot stars in 
the core helium burning stage (hot "horizontal branch" 
or HB stars). These hot HB stars could potentially 
boost the Balmer absorption line strengths sufficiently 
to mim ic younger ages for old globular clusters (e.g., 
iLee fc Worthev 2005: Maraston fc Thomas . 200Qi) . Ide- 
ally, cluster ages should be determined from deep color- 
magnitude diagrams (CMDs) which reveal the main se- 
quence turnoff; to date however, this has only been ac- 
complished for a single cluster beyond the Magellanic 
Clouds (M31-SKHB 312: iBrown e talll200l. using 129 
orbits of Hubble Space Telescope (HST) time. 

In this Letter, we present integrated-light (optical) 
spectroscopy for nine ancient star cluster candidates in 
the nearby spiral galaxy M33. These clusters have deep 
CMDs available from HST/WFPC2 observations, which 
clearly reveal the HB morphology. The combination of 
HB morphology and absorption line measurements result 
in robust relative age estimates, and are ideal to search 
for intermediate-age clusters in M33. 

2. DATA AND ANALYSIS 

2.1. Integrated- Light Spectroscopy 

We obtained integrated spectroscopy for ^150 star 
clusters and cluster candid ates in M33 using the 
MMT/Hectospec instrument (jFabricant et al.ll2005j) at 
a dispersion of 1.2 A pix~^. In addition to clusters, a 
significant number of background regions were chosen, 
in order to sample the range of background light pro- 
vided by the galaxy itself. The clusters were observed 
on the night of October 27/28, 2005, with four 1100 sec- 
ond exposures (for a total of 73.3 minutes). Each of the 
300 Hectospec fibers subtends 1.5" on the sky. The data 
were reduced using the HSRED reduction pipeline devel- 
oped by R. Cool'"^. All observations were bias subtracted, 

^ Available from |http://mizar. as. arizona. edu/~rcool/hsred| 
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Fig. 1. — The blue portion of the background subtracted spec- 
tra of nine M33 globular clusters, taken with MMT/Hectospec are 
shown. The spectra have been normalized by a spline fit to their 
continua. 



brate the index mea surements to the Lick/IDS system 
ijBurstein et aUlToS^ , and instead focus on relative clus- 
ter ages. A comparison of the Une indices measured for 
the nin e clusters pre s ented in thi s work with predictions 
of the Tho mas et all 1120031120041) models is illustrated in 
Figure El 

2.2. Color-Magnitude Diagrams 

The CMDs for the nine clusters analyzed in this work 
were originally presented in Sarajedini et al. (2000; here- 
after SOO). These are based on Vpsssw and /f814W ob- 
servations taken with the HST/WFPC2 instrument. We 
refer the reader to SOO for details of data acquisition, re- 
duction, and analysis. The CMDs show that only M9 
and U77 have blue HB stars; the rest of the M33 clusters 
in the SOO sample have HBs which lie completely red- 
ward of the RR Lyrae instability strip, in the so-called 
"red clump" . 

Here, we present new CMDs based on the SOO data for 
two clusters of particular interest - M9 and CSS. These 
have been radially cleaned, so that blends and other con- 
taminants from the crowded central regions are not in- 
cluded in the CMD. In Figure Owe show the new CMDs 
for these two clusters which are further discussed below. 



overscan corrected, and trimmed. The science exposures 
were flat-fielded and combined together to eliminate cos- 
mic rays, and one dimensional spectra were extracted 
and wavelength calibrated. The final resolution mea- 
sured from the FWHM of HeNeAr comparison spectra 
was 4.7 A. Spectra for the entire M33 cluster sample 
will be presented in a future work. Here, we focus on a 
subset of nine clust ers which have av ailable CMDs reach- 
ing below the HB ( Saraiedini et aLlllQQH. .2000. hereafter 
SOO). 

Because the underlying galaxy light within each object 
fiber can affect the absorption line strengths, we created 
a mean background spectrum for each cluster. This was 
accomplished by averaging several background spectra 
chosen to match the galaxy background level measured 
in an annulus around each cluster from reduced M33 im- 
ages provided by the Local Group Survey^). These cus- 
tomized background spectra were then subtracted from 
each object. For the nine clusters presented here, the 
cluster spectra have such high signal-to-noise, that the 
background subtraction has little impact. In Figure ^ 
we present integrated, background-subtracted spectra of 
the globular clusters discussed in this work. 

The relative ages and abundances of stellar popula- 
tions can be estimated from integrated-light spectra, 
by comparing the absorption line strengths for age and 
metallicity-sensitive lines. The measurements and ac- 
companying uncertainties are computed using the tech- 
niques desc ribed in . Pu zia et al. ( 2002), for li ne indices 
define d in I Wort he vM jl 9941 and Worthey fc Ottavianil 
1)1997(1 . The passband definitions were shifted to ac- 
count for the radial velocit y of each clust e r, ass uming 
the measurements given in 'ChaiLdar ll200l . Be- 

cause no observations of Lick standard stars have been 
made to date with MMT/Hectospec, we do not cali- 

® http:/ /www. Iowell.edu/~masey/lgsurvey. 



3. DfSCUSSfON 
3.1. Evidence for an Intermediate- Age Cluster 

In Figure we show the measurements for the metal- 
sensitive absorption line index [MgFe]' versus the age- 
sensitive Balmer line indices 116 a, Hi a, and H/3, as is 
typically done to estimate the ages and chemical compo- 
sitions of stellar p opulations from integrated-light spec- 
troscopy (e.g., Pu zia et'ani2005ah . For comparison, we 
show simple stellar population (SSP) model predictions 
from Jhomas et al. ( 2004) for six different ages (solid 
lines) and six different values of [Z/H] (dashed lines). We 
also compared the M33 cluster mea surements with the 
solar-s caled SSP model predictions of lBruzual Sz Chariot! 
l|2003j) at 4.7 A resolution (not shown), which yield 
results consistent with the Thomas et al. (2004) mod- 
els. All panels in Figure |21 show a "cloud" of points lo- 
cated towards low metallicities and old ages, although we 
cannot rule out an age spread of the order of At < 0.3, 
which translates into a ~ 4 Gyr spread at 12 Gyr abso- 
lute age. Cluster C38 is clearly offset from the "cloud" 
in every panel, and clearly appears to have a younger age 
than the other clusters - recall that our measurements 
have not been absolutely calibrated to the models, but 
should be robust in a relative sense. Therefore, Figure |21 
implies that C38 is '--^ 5 — 8 Gyr younger than the other 
M33 clusters plotted here. Clusters with Ha seen in emis- 
sion (R14 and U77) are plotted with different symbols in 
FigureEl The Balmer absorption-line strengths for these 
clusters thus represent lower limits. 

Could the stronger Balmer line indices for C38, which 
suggest a younger age, be due to the presence of blue/hot 
HB stars? In Figure |3I we show the resolved CMD of 
C38 based on deep HST/WFPC2 observations for this 
cluster, and also for M9 (one of the older clusters) for 
comparison. C38 shows no blue HB stars, while M9 does. 
Note that the HST observations have sufficient depth and 
resolution to detect blue HBs if they are present, but they 
clearly do not exist in C38. Therefore, the presence of 
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Fig. 2. — Age-mctallicity diagnostic plots for the sample M33 globular clusters, con structed from three different Balmer line indices: H(5a 
{left panel), H7A (middle panel), and H/3 (right panel) vs. [MgFe]' (see also Fig. 9 of lPuzia et al.l2005al . for details). Over-plotted are SSP 
models of Thomas et al. ( 2004) for [a/Fc] = 0.3, metallicities [Z/H]=-2.25, -1.35, -0.33, 0.00, 0.35, and 0.67 dex {dashed lines), and ages 
15, 10, 5, 3, 2, 1 (solid lines). The thick solid line is the 5 Gyr isochrone, and is used to split between old and intermediate age globular 
clusters. Note that some clusters are not plotted in the individual panels due to contamination of their index passband definitions by bad 
pixels. Clusters with Ha observed in emission are indicated by open circles. All clusters are labeled as in Fig. We point out that the 
Lick index measurements are not calibrated to the Lick system so that only relative measures should be derived from the shown diagrams. 



hot HB stars, which has been suggested as a plausible 
explanation for apparent intermediate a ges of clusters i n 
different galaxies (such as in M31; Puz ia et a l] l2005bl) . 
cannot be the correct explanation in this case. The line 
indices and the knowledge of the HB morphology taken 
together, indicate that it is the age of the cluster which 
drives the observed Balmer line strengths. 

There are two features from the resolved CMD of C38 
shown in Figure 13 which add weight to this conclusion. 
First, SOO found that the absolute I band magnitude of 
the red clump for this cluster is more luminous than all of 
the others in their sample, and that this behavior is rea- 
sonably we ll fit by the ^ 2 Gyr theoretical HB models of 
(pi7a7diT999, see also Figure 26 in SOO). Secondly, when 
the CMD for C38 is cleaned to exclude the most crowded 
central regions, there appears to be a hint of a main se- 
quence turnoff (MSTO) in the CMD, as illustrated in 
the left panel of Figure |Sl In contrast, M9 does not show 
evidence for a MSTO. Also in this figure, we show a Z 
= 0.004 isochrone at an age of 2 Gyr from Girardi et al. 
(2000) for comparison. Taken together with the spec- 
troscopic results and the abnormally bright red clump, 
the presence of a MSTO in the CMD of C38 provides 
strong evidence that C38 is a genuine intermediate-age 
(~ 2 — 5 Gyr) star cluster in M33. 

3.2. Properties of M33- 038 

Here we compile known properties of M33-C38. This 
cluster lies at a projected distance of ~ 4.5 kpc from 
the center of M33, has a total T^-band magnitude of 
18.0 1 &nd B - V = 0.73 and V - I = 0.89 col- 
ors (IChri stian fc Schommej 119881 IChandar et alJ 119991 
[2002). SOO estimated that it has an [Fe/H] of -0.65±0.16 
from the slope of the red giant branch, which is con- 
sistent with our spectroscopy given the uncertainties 
(f Fe/H1= -1.1 ± 6). 

IChandar "eTall ll200l measure a radial velocity of 
— 145 ±30 km/s for this cluster. The local Hi (disk) ve- 
locity at the location of C38 is found t o be ^ —100 km/s 
from the radio maps of Warner et a D llT97l . Although 
the cluster velocity measurement deviates from the local 
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Fig. 3.— y-7 versus / CMD for M33 clusters C38 and M9. The 
CMDs have been restricted to include stars in an annulus between 
1 and 2" away from the center of each cluster, in order to minimize 
blending and contamination. A Girardi isochrone for an age of 2 
Gyr and Z = 0.004 is overplotted onto the C38 CMD, which shows 
the red clump stars and a hint of a main sequence turnoff. A 
fiduc ial HB for the Galactic globular cluster M3 I Joh nson & Bolt^ 
119981) is overplotted onto the M33-M9 CMD, which clearly reveals 
the blue horizontal branch stars in this cluser. 



disk motion in Figure 01 of IChandar et all l|200^ . (show- 
ing a comparison of cluster velocities versus the local disk 
motions relative to cluster age, C38 falls in a region which 
is ambiguous as to which M33 component C38 belongs, 
clusters with the most deviant velocities are almost cer- 
tainly in a halo/thick disk component). Therefore, this 
object could belong to either the thin disk or a halo/ thick 
disk component. Implications for these possibilities are 
discussed in the next section. 

Unfortunately M33-C38 does not have any published 
M/L ratios, which would provide a direct estimate of its 
total mass. There are, however, published M/Ly val- 
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ues for three of the clusters presented in this work (H38, 
M9, and E12: lLarsen et alJ l20021. They have an average 
M/Lv = 1.53 ± 0.18, which is consistent with M/Ly 
measurements for Galactic GCs (McLaughlin 2000). The 
integrated spectra presented here corroborate that these 
objects are ancient, even though their exact ages are 
not known. Given the fact that MjLy decreases with 
decreasing age, as an upper limit for C38 we take the 
highest value found empirically bv ILarsen et al.l l|2002D 
of M/Ly = 1.87. As a lowe r limit, we assume the MjLy 
value predicted by the iBruzual fc Chari ot ( 2003) mod- 
els at an age of 2 Gyr and for a low metallicity sys- 
tem (Af/Ly = 1.125). Assuming a foreground reddening 
Eb-v = 0.1 mag towards M33 and a distance modulus of 
24.64 mag, we estimate that C38 has a mass in the range 
5 — 9 X 10^ M0. This is about an order of magnitude 
higher than the masses for old Galactic open clusters, 
and within the range of ancient globular clusters. 

3.3. Implications of Intermediate- Age Clusters in M33 

The confirmed presence of an intermediate-age cluster 
in M33 has important implications for star/cluster for- 
mation processes in spiral galaxies, regardless of which 
structural component C38 belongs to. Below we explore 
the implications of the presence of intermediate-age clus- 
ters in both a thick disk versus halo component, and also 
in the thin disk. 

If C38 belongs to either a thick disk or halo compo- 
nent, it is likely that at least a few of the other clusters 
in our sample belong to the same structural component, 
so that a relatively large age spread among M33 clus- 
ters would be present in that component. In the Galaxy, 
halo GCs are known to be old and have an age spread of 
only ^3 Gyr, whereas the GCs associated with the thick 
disk are essentially coeval ijDe Angeli et al.l[2005[) . We 
suggest that a large age spread for M33 clusters would 
favor a halo over a thick disk origin, for the following 
reasons. The thick disk in the Milky Way is relatively 
metal-rich, and believed to have formed during a single 
event, likely the accretion of a relativel y massive satel- 
lite galaxy w hich puffed up the disk (e.g. iMorrisonllTOQ^ 
lAbadi et all 1200 3: Martin et al. 200|). Such an origin 
for a thick disk in M33 would imply that clusters in M33 
formed for many Gyr relatively undisturbed in the disk, 
until several Gyr ago when an accretion event took place, 
resulting in a thick disk with a cluster population span- 
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ning a large range of ages. Given that to date there is 
no evidence for such an event in M33, and that this type 
of scenario would require that no significant merging oc- 
curred before or since, we suggest that it is unlikely that 
clusters residing in a thick disk would have an age spread 
of several Gyr. 

Therefore, the presence of C38 seems to favor a halo 
component over a thick disk component. This would 
point to a longer timescale for the buildup of the M33 
halo relative to that found for the Milky Way, as sug- 
gested previously by SOO and Chandar et alJ l)2002j) . 

If instead of the halo/thick disk, C38 resides in the 
thin disk of M33 (the velocity measurements are am- 
biguous in discriminating between these two possibili- 
ties), the intermediate age of this object also has im- 
portant implications. First, this would confirm the pres- 
ence of relatively massive, intermediate-age disk clusters, 
which are not widely observed in our own Galaxy (but 
might be hidden behind large columns of dust). Sec- 
ond, the ability of massive clusters (~ 10** Mq) to sur- 
vive in the disk of M33 has b e en su ggested to be lim- 
ited to < 1 Gyr. iLamers et all l)2005|) used cluster sam- 
ples in M33 and their theo ry of cluster dissolution (e.g. 
iBoutloukos fc Lamersl2003l) to estimate directly from the 
cluster age and mass distributions the ability of clusters 
to survive in different environments. For M33, they find 
that 10** Mq clusters will disrupt in log(age) = 8.8 years 
(- 600 Myr), and that a 5 x W Mq cluster wiU dis- 
rupt in ~ 1.5 Gyr. Mass estimates for C38 put it at 
5 — 9 X 10** Mq. Although we do not know its precise 
age, it appears that C38 has likely survived 2 — 3 times 
longer than this prediction. Regardless, if C38 is part of 
M33's thin disk, this points to the formation of a star 
cluster which is more massive than known counterparts 
of a similar age in the Local Group, and also to the abil- 
ity of such an object to survive disruption forces in the 
environment of M33 for several billion years. 

We thank D. Fabricant and N. Caldwell for making the 
wonderful Hectospec instrument available to the commu- 
nity. R.C. and T.H.P. acknowledge support by NASA 
through grants GO-10402-A and GO-10129 from the 
Space Telescope Science Institute, which is operated by 
AURA, Inc., under NASA Contract NAS5-26555. A.S. 
gratefully acknowledges support from NSF CAREER 
grant AST-0094048. 



Goudfrooij, P., Gilmore, D., Whitmore, B. C, & Schweizer, F. 

2004, ApJ, 613, L121 
Johnson, J. A., & Bolte, M. 1998, AJ, 115, 693 
Lamers, H. J. G. L. M., Gieles, M., & Portegies Zwart, S. F. 2005, 

A&A, 429, 173 

Larsen, S. S., Brodie, J. P., Sarajedini, A., & Huchra, J. P. 2002, 

AJ, 124, 2615 
Lee, H.-c, & Worthey, G. 2005, ApJS, 160, 176 
Maraston, C. & Thomas, D. 2000, ApJ, 541, 126 
Martin, N. F., Ibata, R. A., Bellazzini, M., Irwin, M. J., Lewis, 

G. F., & Dehnen, W. 2004, MNRAS, 348, 12 
McLaughlin, D. E. 2000, ApJ, 539, 618 
Morrison, H. L. 1993, AJ, 105, 539 

Puzia, T. H., Saglia, R. P., Kissler-Patig, M., Maraston, C, 
Greggio, L., Renzini, A., & Ortolani, S. 2002, A&A, 395, 45 

Puzia, T. H., Kissler-Patig, M., Thomas, D., Maraston, C, Saglia, 
R. P., Bender, R., Goudfrooij, P., & Hempel, M. 2005a, A&A, 
439, 997 



A genuine intermediate- Age Globular Cluster in M33 



5 



Puzia, T. H., Perrett, K. M., & Bridges, T. J. 2005b, A&A, 434, 
909 

Sarajedini, A., Geisler, D., Harding, P., & Schommer, R. 1998, ApJ, 
508, L37 

Sarajedini, A., Gcislcr, D., Schommer, R., & Harding, P. 2000, AJ, 
120, 2437 

Thomas, D., Maraston, C, & Bender, R. 2003, MNRAS, 339, 897 



Thomas, D., Maraston, C, & Korn, A. 2004, MNRAS, 351, L19 
Warner, P. J., Wright, M. C. H., & Baldwin, J. E. 1973, MNRAS, 
163, 163 

Worthcy, G. 1994, ApJS, 95, 107 

Worthey, G., & Ottaviani D.L. 1997, ApJS, 111, 377 



